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The vortex configurations in the Bose-Einstein condensate of the mixture of two different spin 
states \F — l,m/ = —1 > and |2, 1 > of 87 Rb atoms corresponding to the recent experiments by 
Matthews et. al. (Phys. Rev. Lett. 83, 2498 (1999)) are considered in the framework of the 
Thomas- Fermi approximation as functions of N2/N1, where Ni is the number of atoms in the state 
|1, —1 > and N2 - in the state |2, 1 >. It is shown that for nonrotating condensates the configuration 
with the |1, —1 > fluid forming the shell about the |2, 1 > fluid (configuration "a") has lower energy 
than the opposite configuration (configuration "b") for all values of N2/N1. When the |1, — 1 > 
fluid has net angular momentum and forms an equatorial ring around the resting central condensate 
|2, 1 >, the total energy of the system is higher than the ground energy, but the configuration "a" 
has lower energy than the configuration "b" for all N2/N1. On the other hand, when the |2 > fluid 
has the net angular momentum, for the lowest value of the angular momentum hi (I — 1) there is 
the range of the ratio N2/N1 where the configuration "b" has lower energy than the configuration 
"a". For higher values of the angular momentum the configuration "b" is stable for all values of 
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The realization of Bose-Einstein Condensation (BEC) 
in dilute atomic gases offers new opportunities for study- 
ing quantum degenerate fluids [fij. These condensates 
which contain thousands of atoms confined to microscale 
clouds have similarities to superfluidity and laser, and 
provide new testing ground for many body physics. 

Bose-Einstein condensates are quantum fluids in which 
many particles have the same quantum state. Intrin- 
sic property of the interacting Bose-Einstein condensate 
is superfluidity. Bulk superfluids are distinguished from 
normal fluids by their ability to support dissipationless 
flow. Superfluidity of 4 He atoms has been widely stud- 
ied, however, only recently evidence for critical velocity 
- a key characteristic of superfluids - was observed in a 
Bose condensate of sodium atoms || . 

Superfluidity is closely related to the existence of sta- 
ble quantized vortices which are well known for super- 
fluid 4 He. Vortices can be created in superfluid i He 
by cooling a rotating container with helium through the 
superfluid transition. Quite recently the vortices in one- 
component Bose- Einstein condensate were generated by 
simple "stirring" of 87 Kb atoms condensate ||. Vor- 
tex states in binary mixture of dilute Bose-Einstein con- 
densates have been created 0] using the method pro- 
posed by Williams and Holland |||. Their new scheme 
exploits the possibility of simultaneously trapping oth- 
erwise identical atoms of 87 Kb in two different "hyper- 
fine" spin states |I > and |2 >. (|I > and |2 > denote 
\F = 1, rrif = — 1 > and |2, 1 > spin states of s7 Rb atoms 
respectively). Matthews et al. [|4| applied a microwave 
to the condensate and focused a laser beam at various 
points around its circumference, splitting the atoms into 
two hyperfine states. The result is a double condensate 
with one condensate at rest in the center of the cloud and 
the other in a unit (or multiple) vortex state. Controlling 
the temporal and spatial dependence of the microwave- 
induced conversion of |1 > into |2 > (and vice versa), one 
can directly create a |2 > (or |1 >) state wave function 
having a wide variety of shapes out of a |1 > (or |2 >) 
ground-state wave function [Q or different numbers of 
atoms in each state. 

In |Q two types of vortices have been formed. 
Matthews et al. put the initial condensate into either 
the |1 > or 1 2 > state and then made a vortex in the 
1 2 > or |1 > state. They considered the time evolution 
of the vortices to study their dynamics and stability and 
showed that the dynamics of the 1 1 > state vortices is dif- 
ferent from the dynamics of the |2 > state vortices due 
to different scattering lengths. 

The simultaneously trapped resting condensates con- 
sisting of the s7 Rb atoms in the |1, — 1 > and |2, 1 > spin 
states have been considered experimentally in |6|] . In this 
case the intraspecies and interspecies scattering lengths 
denoted correspondingly as an, eti2, , 0,22 are in the pro- 
portion an : ai2 : 022 = 1-03 : 1 : 0.97 with the average 
of the three being 55(3)A. It has been observed that the 
atoms with the larger scattering length an in the state 
1 1 > form a lower-density shell about the atoms with the 



smaller scattering length 022 [|| ■ In |ij it has been shown 
that when the 1 1 > fluid forms an equatorial ring around 
the resting central fluid this configuration is rather sta- 
ble. Conversely, the |2 > vortex ring sinks in towards the 
trap center and breaks up. 

The aim of this article is to study the relative stabil- 
ity of different vortex configurations as a function of the 
ratio N2/N1, where Ni is the number of atoms in the 
state |1 > and N2 - in the state |2 >. It is shown that 
for nonrotating phase-separated condensates the config- 
uration with the 1 1 > fluid forming the shell about the 
1 2 > fluid (configuration "a") has lower energy than the 
opposite configuration (configuration "b") for all values 
of N2/N1. When the |1 > fluid has net angular mo- 
mentum, the total energy of the system in nonrotating 
container is higher than the ground energy, but configu- 
ration "a" has lower energy than the configuration "b" 
for all N2/N1. On the other hand, when the |2 > fluid 
has the net angular momentum, for the lowest value of 
the angular momentum hi (/ = 1) there is the range of 
the ratio N2/N1 where the configuration "b" has lower 
energy than the configuration "a". For higher values of 
the angular momentum the configuration "b" is stable 
for all values of N2/N1. It should be noted that the ex- 
periment of Matthews et al. Q corresponds to the case 
where configuration "b" is unstable (N2/N1 w 1,1 = I). 
In our calculations we use the parameters corresponding 
to the experiments by Matthews et al. : approximately 
8 x 10 5 atom in a spherically symmetric potential with 
oscillation frequencies of 7.8Hz in the radial and axial di- 
rections for both spin states. Let us describe our results 
in detail. 

The modern theoretical description of dilute BEC was 
originated by the seminal Bogoliubov's 1947 paper where 
he showed that weak repulsive interaction qualitatively 
changes the excitation spectra from quadratic free par- 
ticle form to a linear phonon-likc structure. To describe 
the trapped condensates at T — one can use the Gross- 
Pitaevskii (GP) (nonlinear Schrodinger) equation for the 
condensate wave function This equation appears as 
the generalization of the Bogoliubov theory for the in- 
homogeneous phase. It was widely used to discuss the 
ground state properties and collective excitations in one- 
component BEC @-|n|. 

In order to derive analytic results, some approxima- 
tions must be used. A commonly used one is the Thomas- 
Fermi Approximation (TFA), which ignores the kinetic 
energy terms. It has been shown that in the case of 
one component condensates the TFA results agree well 
with the numerical calculations for large particle num- 
bers, except for a small region near the boundary of the 
condensate P,|Io|. In fact, even for small numbers of par- 
ticles TFA still usually gives qualitatively correct results. 
In some situations the TFA cannot be used to predict 
quantitative features of the binary mixture of BEC. For 
example, the TFA solution considerably underestimate 
the degree of overlap between the condensates, however, 
it provides an excellent starting point of study. 
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Let us first consider the phase separation in binary 
mixture without rotation. 

In the case of a two-species condensate, letting ^(r) 
(i = 1, 2) be the wave function of species i with par- 
ticle number iVj, we can write two coupled nonlinear 
Schrodinger (Gross- Pitaevskii) equations as: 

- -^VVi(r) + \m^ 2 {x 2 + y 2 + A 2 z 2 )^(r) - 
2mi 2 

- Mi^i(r) + Gn|Vi(r)| 2 Vi(r) + Gi 2 |V> 2 (r)| 2 </>i(r) = (1) 

h 2 l 

- 5 V 2 V 2 (r) + -m 2 w 2 (x 2 + y 2 + A 2 z 2 )V> 2 (r) - 

2m 2 2 

- M2^ 2 (r) + G 22 |V; 2 (r)| 2 i/; 2 (r) + G 12 |^i(r)| 2 ^ 2 (r) = (2) 

Equations (|l|) and (^) were obtained by minimization 
of the energy functional of the trapped bosons of masses 
mi and m 2 given by: 



d A r 



2mi 



|Wi(r)| : 



+ -miw?(a: 2 + y 2 + A 2 z 2 )|^i(r)| 2 + 



2m 2 



|V^ 2 (r)p 



1 

+ 2 
G 22 



m 2W 2 (x 2 +y 2 + A 2 z 2 )|V. 2 (r)| 2 
|^(r)| 4 + G 12 |^(r)| 2 |V2(r)| 2 



(3) 



The chemical potentials fii and fi 2 are determined 
by the relations J d 3 r\ipi\ 2 = A^. The trap poten- 
tial is approximated by an effective three-dimensional 
harmonic-oscillator potential well, which is cylindrically 
symmetric about z axis, A being the ratio of angular fre- 
quencies in the axial direction u Z i to that in the trans- 
verse direction A = 0J Z i/uJi. The interaction strengths, 
Gn,G 22 ,Gi 2 are determined by the s-wave scattering 
lengths for binary collisions of like and unlike bosons: 



G, 



Anh au/mi] Gi 



27r?i 2 ai 2 /m, where m 1 



m^ 1 + m 2 x . 

Let us consider now the phase separation due to inter- 
action between two condensates. In our case of |1, — 1 > 
and |2, 1 > we have \m-\ui\ = \m 2 ui\. We simplify the 
equations by using dimensionless variables. Let us de- 

/ r \l/2 

fine the length scale a± = — - — , and define the 

° - 1 - y miwi J 1 

dimensionless variables r = a±r',E — %u)\E' , ipi(r) — 
y/ Nil ' a\tp'i( v ' ) ■ The wave function ip'i(r') i s normalized 
to 1. In terms of these variables the Gross-Pitaevskii 
energy functional takes the form: 



E' = \J dV [ATrlWi 



| 2 +iV 1 ( 2 :' 2 +y' 2 + A 2 z' 2 )|Vy i | 2 + 



+ N 2 (3 2 \V<P' 2 \ 2 + N 2 (x' 2 + y' 2 + AV 2 )|V 2 | 2 + 



-Ar 2M2 eV|^ 2 | 4 - 



2 ^ 2 ^- Nl N^[\ 2 m 2 



a± 



(4) 



Here /3 2 = mi/m 2 = <^ 2 /w 2 and Ui = 8nauNi/a±. Eqs. 
(|l|) and (||) are rewritten as: 



Wi + (x' 2 + y' 2 + aV 2 M 

, 2 , , , 4nai 2 N 2 m 



fi^ip'i +u 1 \i{;' 1 \ tp[ ++- 

a± m 

/3 2 W 2 + (x> 2 + y' 2 + \ 2 z' 2 W 2 - ^ + 



i|^ 2 | 2 ^=0; (5) 



u 2 p 2 \^ 2 \ 2 ^ 2 + + 



4:Trai 2 Ni m 



Vil a ^ = 0; 



(0) 



where fx^ = 2ni/luo\. 

In the TFA, Eqs. @, (g) and (§) can be further sim- 
plified by omitting the kinetic energy. For TFA the phase 
segregated condensates do not overlap, so we can neglect 
the last terms in Eqs. (H), (jE)) and (g), obtaining from 
(||) and (|J), in different regions (to be determined later), 
the simple algebraic equations: 



^(rO^-^-^ + AV 2 )); 

Ui 

l^(rO| 2 = A>(A4-(p' 2 + AV 2 )) 



u 2 ^ 



(7) 
(8) 



Here p' 2 = x L 



y ■ 



From Eqs. (|7|) and (|^) one can see 



that the condensate density has the ellipsoidal form. 

In the case of phase separation, the energy of the sys- 
tem can be written in the form 



E — Ei -\- E 2 



where 



Ei = —huJiNi 



E 2 = -hujiN 2 



\u 2 P 2 



//., - -//_> r i d 3 r'\iP' 2 \ 4 



(9) 

(10) 
(11) 



In order to obtain Eqs. Eqs. (0)-(|) have been 

used. 

To investigate the phase separation in the mixture we 
first suppose that the condensate 1 atoms form an ellip- 
soidal shell about the condensate 2 atoms (configuration 
" a" ) . To determine the position of the boundary between 
the condensates, we use the condition of thermodynamic 
equilibrium fll2]j : the pressures exerted by both conden- 
sates must be equal: 



Pi =P2 



Pressure is given by E 



Pi 



G, 



(12) 



(13) 



Condensate 2 has the form of the ellipsoid with long 
semiaxis q: 



J 2 + X 2 z' 2 



(14) 
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From Eqs. (0)-(|) and one has the equation 

for q: 

(15) 



Mi - q 



Kp 2 



nq 2 , 



where n = sj (aiim 2 )/(a 22 mi). 

Chemical potentials fi[ and p! 2 can be obtained using 
the normalization conditions J d 3 r'|?/4| 2 = J d 3 r'\ifj' 2 \ 2 = 
1 and are given by: 



Mi 



Mi 



M 2 



(l-| 9 ' 3 + f 9 ' 5 ) 2/5 ' 



2/3 2 (^) 5/2 , (Mi) 5/ V 5 



( /i ' 1 )3/2 g /3 

where g = y/jj^q' and 
M° = 



15 



15A«i 
8tt 



2/5 



(16) 



(17) 



(18) 



From equations (|16[)-(|18[) one can determine the chem- 
ical potentials /i'i and p! 2 an( i the semiaxis of the phase 
boundary ellipsoid q as a function of N\ and iV 2 . The 
energy of the configuration "a" E a = E a \ + E a2 is given 
by: 

15 



Sol = 7^1 Ni I Mi 



105 



4 (^)5/2 

o' 3 2 o' 7 
<L_ _ £ /5 + 1_ 

3 5 y 7 



E, 



1 



0.2 



-TlUJ-i No 



M 2 



15 (m'i) 3/2 , 
4 /3 2 (^)5/2 

/2<T _ , , q' 5 , /2 <?' 7 
fey - 2/j 2 Mi — + Mi — 



/3 



(19) 



(20) 



Let us now consider the opposite case when the con- 
densate 1 2 > atoms form an ellipsoidal shell about the 
condensate Jl> atoms (configuration "b"). In this case 
Eqs. (fHf)- ( pc| ) can be rewritten in the form: 

2 / II 1\ 

Qi = MM2 - 9i)> 



Mi 

(m 2 ') 5/2 = 

15 (M 2 Q 3/2 
2 (M?) 5/2 



0)5/2 



H 3 



3/5' 
2^1 



..// „/3 

Mi<h 



M^f 



= 1, 



(21) 
(22) 

(23) 
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105 
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15 { A) 112 y 

4 P 2 (^/ 2 



(24) 



Here /x'i and /x 2 are the chemical potentials in the con- 
figuration "b", qi — ^/Jj^q'x is the long semiaxis of the 
boundary ellipsoid, Ef, is the energy of the configuration 
"b". 

To estimate which configuration is stable, one has to 
compare E a and E . 

To evaluate AE — E a — Eb in general case it is worth 
first to estimate the energy of the phase boundary which 
arises due to gradient terms omitted in TFA. The surface 
energy per unit area, the surface tension, is defined as 
fx = E s /S, where E s is the surface energy, and S is the 
interface area, a may be written in the form 14 lq| : 

1/2 



flUJl 

2V2a 2 



O12 



- 1 



, V a H a 22 

x mm(Nim 2 + n^i 2 ) 1 / 2 



(uiUsNiNi) 1 / 4 x 



(26) 



Taking into account that the surface area of the ellip- 
soid with the semiaxis a±q has the form: 



S = 2na±q z 1 



1 



:10g- 



1 



(27) 



AVA 2 ^! A-VA^T, 

one can estimate the contribution of the surface energy 
E s = aS to the total energy of each configuration. To 
be specific, we will use the parameters corresponding to 
the experiments of Matthews et. al. Q on 87 Rb atoms. 
In this case mi = m 2 , a± = 3.88 x 10~ 4 cm, A = 1, 
N = Ni + N 2 = 0.8 x 10 6 atoms. As mentioned in the 
introduction, an : et 12 : a 22 — 1.03 : 1 : 0.97 with the 
average of the three being 55(3)A 

In Fig. 1 we show the energies of configurations 
"a" and "b" (including the surface energy) E a J[hu>\) 
(solid line) and Eb/(huii) (dashed line) as functions of 
log 10 (N 2 / Ni) . One can see that E a is always lower than 
Ei,. This is consistent with the qualitative assertion and 
experimental observation that it is energetically favorable 
for the atoms with the larger scattering length to form 
a lower-density shell about the atoms with the smaller 
scattering length |ll||l7]]. It should be noted that the 
surface energy is much smaller than the interaction en- 
ergy because the scattering lengths have very close 
values (see Eq. (p6|)). 

Let us now consider the condensates with nonzero net 
momentum. For vortex excitation with angular momen- 
tum Mj , the condensate wave function is given by 

^(r) = |^.(r)|e^. (28) 

After substituting the wave function for the vortex 
excitation (|2|) in Eq. (||), the equation describes the 
ground state of bosons in an effective confinement po- 
tential llh 2 /2m 1 p 2 + l 2 fi 2 /2m 2 p 2 + Vi + V 2 , where Vi = 
m i Lu l (p 2 + \ 2 z 2 )/2 and p 2 = x 2 +y 2 . So within the TFA 
the density of the vortex state has the form: 



K(r')| 5 



1 



mUM-^ + aV 2 ) 



]'■ 
'1 



12 



(25) |^ 2 (r')| 2 = ^ (p' 2 (h)~( P ' 2 + \ 2 z' 2 ) 



p 2 vr 



(29) 



(30) 
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The important new qualitative feature of a vortex in the 
TFA is the appearance of a small hole of radius £, cx 
li/Hi(k), but the remainder of the condensate density 
is essentially unchanged. The fractional change in the 
chemical potentials caused by the vortex (^(h) — 
can be shown to be small [p|]ll|l, of the order of 1/iV 4 / 5 . In 
calculation of physical quantities containing condensate 
density it is sufficient to retain the no- vortex density 
and simply cut off any divergent radial integrals at the 
appropriate core sizes £f = or £§ = /3 2 ^/a*2- Note 

that using the unperturbed density for calculation of the 
vortex properties corresponds to the hydrodynamic limit. 

In the case of the phase segregated condensate, one 
finds from Eqs. (HQ and that the energy 

change due to the presence of the vortices AE(l%, h) has 
the form: 



AE = AE 



iVi 



AE 



1. 



N 2 



2 

2fl2 



i r i 2 r 2 



(31) 



and 



In the hydrodynamic limit ipi is given by Eqs. 
(§• 

Let us consider the configuration " a" . In the hydrody- 
namic limit the location of the phase boundary is given 
by Eq. @. From @ one has: 

AE° Nl _ 5/?(/4) 3 / 2 
(M?) 5/2 



\huj l N l 



1 ~ 3* 



hi 



hi 



i5i 2 M) 1/2 q' 



h 3 

h 

/4 



-|V* 



(32) 



1 



/ 12 . 

ft? I x 



/4<z' 2 



(33) 



The energy for the configuration "b" has the form (the 
location of the phase boundary is given by Eqs. ([2l|-p3|)): 



In- 



2 mu'n 



In 



V 2 <?2 



2(m?)V2 



(34) 



x In 



Ml 



/4Y/ 2 



(35) 



Using Eqs.(|3j^|35|) one can compare the energies of the 
configurations for different values of the vortex excitation 
net angular momenta and different number of state 1 1 > 
and state |2 > atoms. The results of calculations for the 
case of nonzero net momenta are presented on the figures 



2-7. The parameters corresponding to the experiments 
by Matthews et al. 0] were used. 

The figures 2-4 show the situation with nonrotating 
condensate |2 > atoms and a vortex of condensate |1 > 
atoms. In Fig. 2 we show the rotational parts d32]-j35|) of 
the energies of configurations "a" and "b" AE a /{htui) 
(solid line) and AE},/(hu)x) (dashed line) as functions of 
log 10 (iV2/iVi) in the case when the condensate of |1 > 
atoms has net angular momentum I = 1 (lower curves) 
and I — 2 (upper curves). For the configuration "a" 
this means the rotation of the outer shell and for the 
configuration "b" - of the inner part. 

In Fig. 3 the corresponding total energy is shown, and 
Fig. 4 shows the difference between E a and Eb for this 
case. One can see that E a is always lower than Eb. 

Let us consider now the case when the condensate 1 2 > 
atoms has the net angular momentum I = 1 or I = 2. 
This case corresponds to the rotation of the outer shell 
for the configuration " b" and of the interior part of the 
condensate for the configuration " a" . The rotational part 
is shown on the Fig. 5, the total energy - on Fig. 6. The 
difference between the total energies E a and Eb shown 
in Fig. 7 presents the most interesting of our results. 
In the I = 2 case the "b" configuration with vortex ring 
of condensate |2 > atoms in the outer shell is stable in 
the whole range of the relative concentrations Ni and 
N2 while in the case of I = 1 vortex of condensate |2 > 
atoms this configuration is stable only in the regions of 
low N2/N1 or N1/N2 concentrations. When the values 
Ni and N2 approach one another the condensate |2 > 
atoms ring sinks in toward the trap center. 

These results are in agreement with the experiment 
Q at approximately equal values of Nf. when the con- 
densate |1 > has a net angular momentum, it forms an 
equatorial ring around the central condensate |2 >. Con- 
versely, a condensate |2 > vortex forms a ring that tends 
to contract down into the condensate 1 1 > . 

It should be noted that the similar conclusions have 
been made in Ref. fl8|| . Using the dynamic stability 
analysis of the Gross-Pitaevskii equations, the authors 
of Ref. |ll| have investigated the instability mechanism 
of the configuration "b" and concluded that stabilization 
of this configuration can be attained by the controlling 
the relative population of both species. 

To summarize, we have shown that for nonrotating 
phase-separated condensates the configuration with the 
|1 > fluid forming the shell about the |2 > fluid (config- 
uration " a" ) has lower energy than the opposite configu- 
ration (configuration "b") for all values of N2/N1. When 
the |1 > fluid has net angular momentum, the total en- 
ergy of the system is higher than the ground energy, but 
the configuration " a" has lower energy than the configu- 
ration "b" for all N2/N1. On the other hand, when the 
1 2 > fluid has the net angular momentum, for the low- 
est value of the angular momentum hi (I = 1) there is 
the range of the ratio N 2 /Ni where the configuration "b" 
has lower energy than the configuration " a" . For higher 
values of the angular momentum the configuration "b" 
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is stable for all values of N2/N1. 
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FIG. 3. The total energies of configurations "a" and "b" 
E a /(Tiu)\) (solid line) and Eh/{%ui) (dashed line) as func- 
tions of log 10 (A r 2/A/"i) in the case when the condensate of |1 > 
atoms has net angular momentum I = 1 (lower curves) and 
/ = 2 (upper curves). 

FIG. 4. The difference of total energies for configurations 
"a" and "b" (E a ~ E b )/(huji) as a function of log 10 (7V 2 /iVi) in 
the case when the condensate of |1 > atoms has net angular 
momentum I — 1 (solid line) and I = 2 (dashed line). 
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FIG. 5. The rotational energies of configurations "a" and 
"b" AE a /{Tiu>\) (solid line) and AEb/{huji) (dashed line) as 
functions of log 10 (A r 2/A r i) in the case when the condensate of 
1 2 > atoms has net angular momentum I — 1 (lower curves) 
and I = 2 (upper curves). 



FIG. 6. The total energies of configurations "a" and "b" 
E a /(huii) (solid line) and Eb/(TiuJi) (dashed line) as func- 
tions of log 10 (-/V2/iVi) in the case when the condensate of |2 > 
atoms has net angular momentum I = 1 (lower curves) and 
I = 2 (upper curves). 



FIG. 7. The difference of total energies for configurations 
"a" and"b" (E a -E b )/(htJi) as a function of log 10 (7V 2 /iVi) in 
the case when the condensate of |2 > atoms has net angular 
momentum I — 1 (solid line) and I = 2 (dashed line). 



mat/9912308; Phys. Rev. Lett. 84, 4264 (2000). 



FIG. 1. The total energies of configurations "a" and "b" 
E a /(huiiN) (solid line) and Eb/(huiN) (dashed line) for non 
rotating condensate as functions of log 10 (./V2/iVi). 



FIG. 2. The rotational energies of configurations "a" and 
"b" AE a /(huJi) (solid line) and AEb/(hcui) (dashed line) as 
functions of log 10 (iV2/iVi) in the case when the condensate of 
|1 > atoms has net angular momentum / = 1 (lower curves) 
and I = 2 (upper curves). 
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